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Abstract 

For a subquadratic symbol H on x — r*(R''), the quantum propagator of the time 
dependent Schrodinger equation = Hip is a Semiclassical Fourier-Integral Operator when 
H = H{x,hDx) (ft-Weyl quantization of H). Its Schwartz kernel is describe by a quadratic 
phase and an amplitude. At every time t, when h is small, it is "essentially supported" in a 
neighborhood of the graph of the classical flow generated by H, with a full uniform asymptotic 
expansion in h for the amplitude. 

In this paper our goal is to revisit this well known and fondamental result with emphasis on 
the flexibility for the choice of a quadratic complex phase function and on global estimates 
when h is small and time t is large. One of the simplest choice of the phase is known in chemical 
physics as Herman-Kluk formula. Moreover we prove that the semiclassical expansion for the 
propagator is valid for \t\ << ^| log/i| where 5 > is a stability parameter for the classical 
system. 



1 Introduction and Results 

Let us consider the time-dependent Schrodinger equation 

^^d^^g 'A(^ = io)=^o, (1.1) 

ot 

where ip is an initial state, H{t) is a quantum Hamiltonian defined as a continuous family of self- 
adjoint operators in the Hilbert space L'^{W^), depending on time t and on the Planck constant 
h> which plays the role of a small parameter in the system of units considered in this paper. 
H{t) is supposed to be the ?i-Weyl-quantization of a classical smooth observable H(t,X), X = 
(a;,^) G R"* X R'* (see [27] for more details concerning semiclassical Weyl quantization). 
Our main results concern subquadratic hamiltonians H; that means here that H{t, X) is continuous 
in teR,C°° smooth in X G M^d ^^^^ satisfies, for every 7 e N'^'^, \j\ > 2, 

\dlH{t,X\<CT,-y, '^t,\t-to\<T, VXeM^d (^2) 
where dx = and Ct.j > 0. 

Let us introduce some classes of symbols ("classical observables" ) defined as follows. Let be 
m,n G N. 

Definition 1.1 We say that a symbol s is in Om{n) if s is a smooth function on the Euclidian 
space K" such that for every 7 G N", I7I > m we have 

|sU,7 := sup \dlsiX)\ < +00 (1.3) 

xeR" 
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If s{e) depends on a parameter e €z P we say that s{e) is bounded in Orn{n) if for every 7, we have 

sup |s(£)|oo.7 < +00. 

eeP 

It is well known that the subquadratic assumption entails that equation (jl.ip is solved by a 
unique quantum unitary propagator in L^(M'*) such that ipt — U{t,to)ipo, Vt G M. For the same 
reason, the classical dynamics is also well defined Vt E R. zt — {qt,Pt) is the classical path in the 
phase space M.'^'^ such that Ztg = z and satisfying 

( 4t= dpH{t,qt,pt) ^ 
\ Pt = -dHq{t,qt,pt), qto q, Pu, ^ P 

It defines an Hamiltonian flow : (t>^{z) ~ zt = z). Let us introduce the stability Jacobi 

matrix of this Hamiltonian flow : F(t) — dz(t>*{z). F{t) is a 2(i x 2d symplectic matrix with four 

dy. d blocks, F{t) = ( ^* ^* ) , where 
y-'t F>t J 

. dqt dqt dpt dpt , . 

At = -5-, Bt = 6i = — , Dt = -TT- (1-5) 
oq op Oq op 

We also introduce the classical action 

Sit,z)= f {p,-qs-H{s,Zs))ds (1.6) 

J to 

where u ■ v denote the usual scalar product for u,v £ K'^, and the phase function 



$(t, Z] X, y) = S{t, z)+pfix-qt)-p-iy-q) + ^{\x - qt\^ + \y - q\^) 



(1.7) 



For applications it is useful to introduce semi-classical subquadratic symbols. These symbols have 
an asymptotic expansion in the semiclassical parameter h > 0, H^{t,X) x ''^^h-' Hj{t, X) such 

that the following conditions are satisfied. 

Vj > 0, Hj{t, •) e C(2-i)+ (2d) and are bounded in 0(2_j)+ (2d) for t E R, (1.8) 
\fN>l,h-^-^{H{t,X)- ft^iJj(t,X)) is bounded in Co for te Hand ?ie]0,l]. (1.9) 

0<j<Ar 

Let us recall the definition of Weyl quantization. For any symbol s in 0„i{2d),a,nd for any ip E 
i"^), we have 

OpjrW^(x) - i27rh)-'' f f e*(--^)-«s(^,^)^(y)dyde (1-10) 

././B2d Z / 



x + y 

We shall also use the notation s — Op^ [s] . 

The Herman- Kluk formula is included in the following asymptotic result which will be discussed 
in details in this paper. This formula was discovered by several authors in the chemical-physics 
litterature in the eighties. We refer to the introductions of [12] and [2£| for interesting historical 
expositions. It is rather surprising that until the recent paper [29] there was no explicite connexion 
in the mathematical litterature between the Herman-Kluk formula and Fourier-Integral Operators 
with complex phases. 
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Theorem 1.2 Let be H^{t) a time dependent semiclassical subquadratic Hamiltonian and K^{t\ x, y) 
be the Schwartz kernel of its propagator U^{t, to). Then there exists a semi- classical symbol of order 
0, a^{t;z) = aj{t;z)h-' where Oj is continuous in t, 

a<j<+oo 

K^{t;x,y)^ I eTi'^^'^''-^'-''y^a{h;t;z)dz (1.11) 
in the LP' uniform norm. More precisely, if we denote 

K^''^^\t-x,y)^{2nh)-'^'^'^ ( e^^^'-'''^^y\y^ a,{t-z)h')dz (1.12) 

J^^' 0<j<N 

and U^'^''^\t,to) the operator, in L^(M''), with the Schwartz kernel K^^'^^t^Xjy), then, for every 
T > and every N > 1, there exists C{T, N) > such that for the L^ operator norm we have 

\\uh{t,to) - u^'''''\t,to)\\ < c{T,N)n''+\ yt,\t-to\<T, ne]o,i]- (i.i3) 

The leading term is 

ao{t;z) = det^^\At+Dt + i{Bt-Ct))exp(-i f Hi{z,)ds) (1.14) 

Jta 

where the square root is defined by continuity starting from t — to (ao{to; z) = 2'^^'^). 

Moreover, the amplitudes Oj are smooth functions defined by transport equations (see the proof 

below) and, for every T > they are bounded in Oq for \t\ < T . 

In ^9J the authors give a rigorous proof of this result with an additionnal hypothesis : they 
assume that _ff(x,^) is a polynomial in ^. Here we consider more general subquadratic symbols. 
In particular our result applies to relativistic Hamiltonians like ■\/l + + V{x). Using a global 
diagonalization (see [28], section 3), the result can be extended to Dirac systems. 
Similar results are true with more general quadratic phases and for systems with diagonalisable 
leading symbols (see [HHQ). Let us define the quadratic phase 

^^'^^^^\t,z-x,y) = S{t,z)+pf{x~qt)-p-{y-q) + ]^{Qt[x~qt)-{x-qt)-f{y-q).{y~q)) (1.15) 

where F, Qt are complex symmetrix matrices with a definite-positive imaginary part, Qt is in 
t. r is constant, O* may depend smoothly on t and z such that the following condition is satisfied 

3cT > 0, '^Otv.v > —\vP, Vt, lil < r, Vz e R^"^ (1.16) 

CT 

V7,|7| > 1,3Ct,^, \\d2et\\ < CT.rn e R"" Mt\ < T. (1.17) 

So we have 

Theorem 1.3 Under the assumptions of Theorem I j.^l and iL16\) . Ijl.lT^ , we have 

K{t-x,y)^(2Tihy^'^'^ [ e^**^"^'(*'"'"'^)/(a;t;z)dz (1.18) 



where f{h;t;z) = fj{t\z)h^ with the same meaning as in Theorem II .21 

0<j< + cxj 

In particular 

fo{t, z) = 2'*/2deti/2[M(et, r)] (1.19) 

where 

M(et,r) = i{c + DT-e{A + BT). 
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There exist several methods to prove this theorem. In |29j the authors prove it as a consequence 
of a symboHc calculus for FIO with complex quadratic phases. In |5] the authors proved a weaker 
result for F = il and Qt = Ft is determined by the propagation of Gaussian coherent states: 
Ft = {C + Dr){A + Br)~^ (see section 2 of this paper). Laptev-Sigal in |^ have also considered a 
similar formula for the propagator (see section 5 of this paper) but assume that the initial data has 
a compact support in momenta. Kay [22| explains how to compute all the semiclassical corrections 
ttj but did not give estimates on the error term, so its expansion is not rigorously established. 
Here we choose another approach, may be more explicit and simpler. We shall prove the general 
theorem 11.31 as a consequence of the particular case of Theorem 11.21 by using a real deformation 
of the phase on the simpler one Moreover we give a direct proof of Theorem 11.21 

proving the necessary properties for Fourier integrals with complex quadratic phases. This way we 
can get easily explicit estimates for the error terms for large times. 

Let us assume that conditions on H{t) are satisfied for T = +oo. Moreover assume that there 
exists a positive real function fJ,{T) > 1, T > 0, such that the classical flow 0* satisfies, for every 
multiindex 7, I7I > 1, we have for some > 0, 

|9>*'*V)| < C^m(7')'''', for \t\ + \t'\<T, VzeM^'^ (1.20) 

We have discussed in [5] the condition (|1.20p . In particular this condition is fulfilled with /x(T) = 
e'^^for5= sup \\Jd]i xH{t,X)\\. 

Theorem 1.4 Choosing the phase as in theorem \1.2l for j > the amplitudes aj{t,z) satisfy the 
following estimates, for every multiindex 7 there exist a constant Cj^^ such that 

\d2aj{t,z)\ < Cj^\det^/^Mt\fJ.{t)^'+^-^^, G M, Vz e M^''. (1.21) 

Hence we have the following Ehrenfest type estimate. For every N > 1 and every e > there exists 
Cn.e such that we have 

|lC/(t,to)-f/(^)(i,to)ll <CAr,,r(^+i), yt,\t\>^,Whe]0,l]. (1.22) 

In previous works an Ehrenfest time Te — c\ogh~^, c > 0, was estimated for propagation of 
Gaussians in [Hj and propagation of observables in [5] . For Gaussians we got c~ ^ , for observables 
c = In I29j the authors gave an Ehrenfest time without explicit estimate on c. 



2 Gaussians Coherent States and Quadratic Hamiltonians 

The phase functions c&t^''"' in p.7p and (jl.lSp are closely related with Gaussian coherent states. 
This can be seen by proving a particular case of Theorem ) 1.21 for quadratic time-dependent Hamil- 
tonians: 

Ht{<l,P) = ^{Gtq-q + 2Ltq-p + Ktp-p) 

where q,p & K'^, Kt, Lt, Gt are real, dxd matrices, continuous in time i G K, Gt, Kt are symmetric. 
The classical motion in the phase space is given by the linear differential equation 

where L'^ is the transposed matrix of L, J defines the symplectic form a{X,X') := JX ■ X', 



5 



This equation defines a linear symplectic transformation, Ft, such that -Fb = I (we take here 
to — 0). It can be represented as a 2(i x 2d matrix which can be written as four d x d blocks : 

« = ( « n. ) ■ f^-^^) 

The quantum evolution for the Hamiltonian H{t) is denoted by U{t) (U{Q) = I). We can compute 
the matrix elements of U(t) on the coherent states basis tpz- This has been done in Littlejohn 
[24] (p. 249, (6.36)), Bargmann 3 , Fedosov [12], [10]. We follow here the presentation given in 
[10] . Let us introduce some notations which will be used later, g denotes the Gaussian function: 
g{x) = 7r~''/'*e~l^l and is the dilation operator Aritp{x) = h^''-/'^tp{h^^/^x). So tpo ~ A^g, and 
the general Gaussian coherent states are defined as follows. 

^P=f{z)^^^\ (2.25) 

where T{z) is the Weyl translation operator, z = {q,p), 

f{z) ^ exp (^^{p ■ X - q ■ hD,)^ (2.26) 

where = —i-§^ and z = {q,p) G M'^ x W^. ip^^'> is the Gaussian state: 

(^(r) (a;) ^ (7rft)-''/4ar exp (^l^Tx.x^ (2.27) 

where F is a complex symmetric matrix such that 5F is definite-positive, ar is a normalization 
constant, (ar = det^^^^F). 

It is convenient to introduce here the Siegel space S+(c?) of c? x c? complex matrices F such that 
SF is definite-positive, (see in [13] properties of S+(c?)). 

Let us define the Fourier-Bargmann transform t'"^^ as follows, "0 G L^(M''), 

4^)[^](z)^(27rft)-^/2(0,^(r)). (2.28) 

the following coherent state living at z, z — {q,p) e M'* X M'*, X e R"^, 

)(.) . (.;.)-/^.r exp(l(p • X ~ ^) + !£(fL^|_(:I^) (2.29) 

JFg"'' is an isometry from L^(M'^) into i^(R^'') (with the Lebesgue measures). If F = il we denote 
J^e = -^g; its range consists of F e L^(]R.^'^) such that exp -I- F{q,p) is holomorphic in 
C'' in the variable q — ip. In other words, 

J^B^z) = E^iq ~ zp)exp (^-^ - (2.30) 

where E^, is entire in (see [25]). Moreover we have the inversion formula 

^[x)^ \ J^^^^[i/'](z)^f )(x)dz, in the - sense. (2.31) 



and — ip. 



These properties are well known (see [251 [5]). Sometimes we shall use the shorter notation tp^ 
r-(r 

B 
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Let us denote by i?[^t] the quantum propagator for the HamUtonian H(t) (this is the metaplec- 
tic representation of Ft) and iiT^^*) its Schwartz kernel. We know that AfiR[Ft]g is the following 
Gaussian state [TUl [H] , 

AnR[Ft]g{x) = {Trh)-''/W{t) cxp (^-^^tx.x^ (2.32) 

where ar{t) = [det(yli + TBt)]^^^'^ar, the complex square root is computed by contimiityQ from 
t = to ^ 0, and 

Tt = iCt+TDt)iAt+TBt)-\ Tt,^T. (2.33) 



Proposition 2.1 We have the following exact formula 

/f(^')(x,t/)-2'^/2(2^n)-3'^/2deti/2fM(^^) / e*<^'''(*'^^-'^)dz (2.34) 

V I ) J^2d 

where T, St G ^+{d), Qt is in t; Af(ef , T) = C + L»f - Of (A + BY) and 

<^'^^*-^'>{t,z;x,y) = ^{qfPt-q-p)+Pf{x-qt)~p-iy-q) + ^{'dt{x~qt)-{x~qt)-r{y-q).{y-q)) 
Let us remark that here the action is S{t, z) — ^{qt ■ pt — q ■ p) ■ 

First of all let us remark that the integral (I2.34p is an oscillating integral and is defined, as usual, 

by integrations by parts. We shall give two proofs of this formula. 

Proof I. We start with any Fq in the Siegel space E+ [d) . Using the formula 

i^{x)^{2Tih)-'' ( (^,^^o^^ro^^ 

we get the formula 

K(^^\x,y) = {27rhr'' f ^^F" {y)^f:\x)dz (2.35) 

So, we get 

if(^*)(x,y) = (27rft)-3'^/2^o(i) / e^*'"""°'(*'^'^'^)dz, (2.36) 

where 

dei^/^{A + BT^) 
Now we shall transform the phase $('"'^'"") into the phase $(®''"o). 

Let us introduce Q{s) = sG + (1 - s)Ft, < s < 1. We have e(s) G ^+{d). We want to find k{t, s) 
such that A:(<, 0) — ko{t) and 



d_ 

We have 



{k{t,s) I e^*'^"^"'(*'^'"'^)dz) =0, Vse [0,1]. (2.37) 

^ ./Bad ^ 



^e«* =_(et-Ft)(a;-gt)-(a;-'zOe« 



* this definition of det^'''^ is difl'erent that the det^'''^ function on E+(d), this is explained in |10| to compute 
Maslov index 
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The main trick used here and later in this paper, and also in all the previous papers on this subject 
([231 [2H US])) is to integrate by parts to convert each factor {x — qt) into H, using the following 
equality 

{dg + f <9p)$®^^ = (C" + TD^ ~ [A^ + fB^)e) (x - qt) (2.38) 
where denotes the transposed matrix of A. Let us introduce the matrix 

M = M{e, r) = c + Dr- &{a + bt) 

So we have 

Ar(x-gt)et*'^''' = !^{dy + Tdp)eT'^^'". (2.39) 

Let us remark that M is invertible. This is a consequence of the following Lemma (see [11], [13] 
or ^5], appendix A, for proofs). 

Lemma 2.2 For every linear symplectic map in F : T*(R'^) T*(R'^), F = and every 

r e i;+ (d) , {A + by) , (C + DT) are non mvertible in and (C + DT){A + BT)-^ e T.+ (d) . 
So we have 

M = C + L»r - e{A + BT) ^ {{C + DT){A + Bvy^ - e) (A + Bvy^ 

But (C + DT){A + BT)-^ - 9) G S+(d) so is invertible. 
Denote M(t,s) = M{Qs,^t)- Let us recall the Liouvillc formula 

dsdet{M{t,s)) ^ det{M{t, s))TT(dsM{t, s)M{t, s)-^y (2.40) 

So, integrating by parts in (g, p) we get 

, / X , / „xdet^/^M(i,s) , , 

Now we have to compute ^^^.1/2'^^^ ■ A simple computation gives M{t, 0) — {D — rtB){To — Fq). 
The proof of (|2.34p follows from the formula 

det(i:>-r(S) =det(A + Bro)"^ (2.42) 

This equality follows from the symplecticity of F {D^ B = B^D). We have B^'TtB - B = 
— {A + Bro)^^B. So we get (|2.42p if deti? ^ 0. The general case follows by a density argument. 
Let us remark that can exchange the role of Q and F by considering the adjoint U(t)* of U{t). □ 
Proof II 

We solve directly the Schrodinger equation 

^^n^^-H{t))ij{t,x)^0 (2.43) 
for any initial data ip{x) := ip{0,x), ip S 5(R'') using the ansatz 

i>{t,x)^{2TThy^'^/^k{t) j e**'''''''(*'^'^'^V(2/)rf2dy (2.44) 



We have to compute k{t) such that fc(0) = 2^^/^. Let us remark that if we integrate first in y then 
the integral (|2.44|) in z converges because the Fourier-Bargmann transform of '0, ^B^i is in the 
Schwartz space S{M.'^'^). 

For simplicity we assume here that = F = il. The general case can be reached by the same 
method or by using the deformation argument of proof I as we shall see later for more general 
Hamiltonians. 

Here the Hamiltonian H(t) is a quadratic form. So using dilations we can assume that h = 1. A 
simple computation left to the reader, gives the following 



8 



Lemma 2.3 

{g-^H{t)g){x) =Gx-x + i{L + U)x ■ x - Kx ■ x + Tr(A' - iL) (2.45) 

where g{x) = e~4~. 
So we get 

(idt-H{t))^{t)^{2Tih)-'^'^/'' If e^**^'^'(*'^;"'^)6(t,x,z)V'(2/)dzdy (2.46) 

where 

h{t, z, z) = idtk{t) - k(t){E{x - qt) ■ {x - qt) + Tt{K - iL)) 
As in proof I, we integrate by parts in the variable z e M^'', using 

{dg-tdp)<i> = Arix-qt) 

with AI = C — B — i{A + D), which is invertible (see below Lemma [3.2^ . Using the Hamilton 
equation of motion we get 

M = -E{A - iB) - i{K - iL)M. (2.47) 
So, we find the following differential equation for k{t), 

k^^TT{{MM)k. (2.48) 

Using the Liouville formula, we get again (|2.34p for this particular phase. □ 

3 Proof of Theorem 11.21 and Theorem 11.41 

As usual for this kind of problems there are two steps : 1-Determine the amplitudes aj solv- 
ing by induction transport differential equations, 2-Estimate the error between the approximated 
propagator and the exact one. 

3.1 Transport equations 

It is convenient to write 

ei* = {TTtif/^^,,ix)ip,{y)ei'^^'^''''>+'^P-''-P'-^''>^^^ (3.49) 

Then we have to compute H^{t)(pzt ■ It is not difhcult to add contributions of the lower order terms 
of the Hamiltonian, so we shall assume for simplicity that H^{t) = HQ{t) := H{t). 

Lemma 3.1 For every N > 2 we have 

m^zAx) = ^dlH{t,zt)nJ^^)^,,{x)+fi^''^'^^'T{zt)AnOpnRN{t,ztMx) 
h\<N ^' ^ 

(3.50) 

where 

RNit,Zt,X)^ / ^ dlH{t,zt + sVhX)XUs (3.51) 

and is a universal polynomial of degree < I7I which is even or odd according \ j\ is even or odd. 
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Proof. Let us recall that = T(z)Ang. In this proof we put zt ~ z. An easy property of Weyl 
quantization gives 

tq^f{z)U{t)f{z)kn = Op'^iHiVh^+z)] (3.52) 

So the Lemma follows easily from the Taylor formula with integral remainder. □ 

In this first step we don't take care of remainder estimates, this will be done in the next step. 
Let us denote I(a, $) the formal operator having the Schwartz kernel 

Kaix,y)^{2TThy^'^/^ f e^*(*'^'"'f)a(t,z)dz (3.53) 

From the Lemma 13.11 we can write 

H{t)I{a, $) - I{b, $), where 

^ - E ^9lH{t, z,)n, (^) a (3.54) 

We have 



/3<7 

The quadratic part can be computed as for quadratic Hamiltonians and the linear part disappears 
with the classical motion. So we have 

b - H{t, zt)a + {dqH{t, Zt) + idpH{t, Zt)) ■ (x - qt)a + 

n(i?(^) . (^) + Tr(A' - ^L))a (3.56) 

where we denote dx x^i^^ the Hessian matrix of H{t). We have 

dj,^xH{t,zt)^ E^G + 2zL-K. (3.57) 



with G := dlgH(t, Zt), L := dl^H{t, zt), K dl^H{t, zt) 

Here the stability matrix Ft = ^* j satisfies Ft — Jdx x-^i^^ ^t)Ft, Ft=Q = I. 



I nauiniicn ±-t — J i^^ 

As in the quadratic case we want to transform the power of {x — qt) into power of h. 



Lemma 3.2 Let us denote Mt — {Ct — Bt) — i{At + Dt). We have 

\AeiMt \ > 2-^ and (3.58) 

Hdq - iap)e** = iM[{x - qt)eT^'^ (3.59) 

Proof. For simplicity, let us forget the lower index t. 
Let us consider the 2d x 2d matrix 

l + F + ^J(l-F)-(^ + ^-'^' B + ^iI-D)\^I + A-^C -^iD + ^B) + ^\ 
Using the Lemma 4 in [13 , Appendix A, we get 

det{I+F+iJ{I-F)) = dct{{l+A-iC){l+D+iB)-{A-iC-I){D+iB-I)) = 2"^ dot {A+D+i{B-C)) 

(3.61) 

Using that F is symplectic, we get 

{1 + F + tJ{I - F))* (I + F + tJ{I - F)) = (I + F^)(I + F) + (1 - F^)(I - F) > hd (3.62) 
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hence (|3.58p follows. 

Let us recall classical computations for the derivatives of the action 

d,S = {dqqtYpt-p (3.63) 
dpS = (dpqtYpt (3.64) 

Then we can compute 9^$, 9^$ and we get p.59|) . □ 
Integrate by parts like in the quadratic case, we get 

{ihdt - H{t))I{a, $) X $) (3.65) 

where 

/ - ih{dta - iTr(AfM-i)a 



+ E ^dlH{t,z,)n,[^)a (3.66) 

|7|>3 ^' ^ 

Hence using the Liouvillc formula, we get the first term 

ao(t, z) = 2'^/2^et^/2(iM) (3.67) 

We shall obtain the next terms Oj by successive integrations by parts. This is solved more explicitly 
with the following Lemma. 

Lemma 3.3 For any symbol b G 00(2^), and every multiindex a G N^'' we have 

[ {x-qtTei'^b{z)dz= V h\^\ j f^^!3{t,z)ei^d^b{z)dz (3.68) 

^<l/3|<l"l 

where fa.pit^z) are symbols of order 0, uniformly bounded in Oo(2d) on bounded time intervals. 
They only depend on the classical flow (j)*{z) and its derivatives. 

More precisely, let us assume that there exists a non positive function fJ,(T) such that for every 
7 G N^*^ we have 

sup \d](t>\z)\ < C^m(7')''^' (3.69) 

\a\<T 

Then we have 

|a|/„./3(z)| <a,;3;,A^(r)i"i-i^i+i^i (3.70) 

Proof. The Lemma is easily obtained by induction on |a| using Lemma [3.21 

Now, to determine the transport equation, we solve inductively on j > 0, the equation 

{ihdt - H{t))l(^ = 0{W+^) (3.71) 

0<fe<j+l 

Reasoning by induction on j > 0, we get the transport equation for aj_|_i(t) by cancellation of the 
coefficient of h^+^ in (|X7T|l . 

dtaj+i{t, z) = ^Ti{MM-^)aj+i{t, z) + bj{t, z), aj+i(0, z) ^ 0, (3.72) 

where 

bj{t,z)= F,^kAt,^)d>k{t,z). (3.73) 

|a|+2fc<2(i+2) 
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Moreover, Fj^k,a{t, z) depends only on the classical flow 4)^{z) and its derivatives and satisfies 

\d:F,^uAt, ■A\ < Q,;c,a,^/i(T)'^-'"'+'^+'^'"'"' (3.74) 
where Cj^k,a,-f only depends on sup \H{t)\oo,-y, 2 < I7I < j + 2. 

|t|<T 

So we get, for every j > 0, 

aj+i{t,z)^ [ det^/^ {M{t,z)M{s,z)-^)bj{s,z)ds (3.75) 



Moreover, from (|X75|) and ^1^, we get the following estimate, for every j > 0, \t\ <T, z e M^'', 

\d]aj{t,z)\ < Cj^^\det^^^Mit, z)\fi{T)*^+\'^\ (3.76) 
with the same remark as in p.74p for the constant Cj^j. 

3.2 Error estimates 

Let us denote 

7^^^(^) = {ihdt - ff(t))l(a(^)(t), (3.77) 

where a'^^\t) = h'^ak- Using Duhamel formula we have 

0<fc<Ar 



\\U^(t)~U"'''{t)\\ <h-' \\n{s)\\ds (3.78) 
where to = 0, U'^it) = U'^it, 0), C/^''^(i) = x(a(^)(i), . 

So we have to estimate ||7?,Ar(t)||. Let us denote K''^\x^y) the Schwartz kernel of Ti-N{t) and 
K''^\X^Y) the Schwartz kernel of 7iN{t) in the Fourier-Barmann representation : 

K^'''>{X,Y)^ ff K^^\x,y)^x{y)Wi^dxdy. (3.79) 



Let be TZN{t) the operator with Schwartz kernel K^^\X, Y). The following Lemma is well known. 
Here we forget N and t for simplicity. 

Lemma 3.4 We have the l? norm estimate 

||7^||L2(R.) < (27rri)-^|l^|U.(M.) (3.80) 

In •particular we have 

||7^||i2(Rd) < (27rn)-'*max|sup J \K{X,Y)\dX,snp J \K{X,Y)\dY^ (3.81) 

Proof For inequality ((3^ we use that the Fourier-Bargmann transform is an isometry. 
Inequality ()3.8ip is known as Carleman (or Schur) estimate. □ 
Using Lemma l3. II we get 



ifW(X,y) = 2-3^/2(7r?i)-'' / {fizt)At,Opr[RN{t)]g,^Y){^x,^M''Ht,z)ei'^''^dz (3.82) 
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where 5{t, z) = 3(1, z) + EliM. 

Using Weyl commutation formula we have 



(^x,^.> = exp(-^^ + ^a(X,z)) (3.83) 
{T{zt)t.nOv'"A^Nit)\g,^Y) = {Opi' [Rn {t)]9 , g r^) . (3.84) 



2d 



We know the Wigncr function Wa.z of the pair {g,gz), Z e R^'' [28] 

W^o.z(^) = 22'^exp(-|X - ||2 - ia{X, Z)) (3.85) 
By a well known property of Weyl quantization [13] . for any symbol s, we have 

{Opr[s]g,gz) = i2nr' [ siX)Wo,z{X)dX (3.86) 



We shall use the following Lemma 

-"1 



Lemma 3.5 Let be f G Oa{2d). For every 7 e N^'' and m > i/iere exists C-y,m smc/i t/iat 



/ X^/(X)e-l^-^l'-^-^^-^dX <C^,™(1 + |Z|)-'" sup \d^fiY)\ (3.87) 

JR2d |Q|<m+|7|; YeK^'' 

Proof It is enough to assume \Z\ > 1. We integrate m times by parts with the differential operator 

2{X^Z)-UZ^ 
^ = 4|X-zP + |JZ|2 (^-^^^ 

using that (r^)™ = ^ lm,o.d'^, with < C„,„(|Z| + \X - Z|) — , where 

I a I < 

6'(X) ^-\X - Z\^ - iJZ ■ X. □ 

So using Lemma [331 we get the following estimate : for every A''; N' there exists Cn.n' (depending 
only on semi-norms |iL(t)|oo,7, 2 < I7I < TV + TV', such that for X, F e M^'^ and |t| < T we have 

|xW(X,y)|<CA.,A..(M(r))^+^'ft^-'^ / e-^(l + ^l^)"'^'|aW(t,z)Mz. (3.89) 

Let us denote 0** = <^°'* = (0*)~^. We have the Lipchitz estimate, for |t| < T, 

W^'Y -z\<ii{T)\Y -zt\ (3.90) 

So we get 



Y -zt\ 



1 + ^ dz 



-N' 



<Cjv' 1 + ^^— ^ (3.91) 



-N' 



and 



|i^(^HX,F)| <C^,^,(MT))^+^'n^(l + ^^^^J)"'" sup |aW(t,z)| (3.92) 



2eR2'i,|t|<T 



Then using Lemma 13.41 and choosing N' > 2d, we get the following uniform estimate for the 
remainder term, for \t\ < T, 

l!7^^^(^)|| < Cjv(/^(r))^+in(^+i)/2 sup |aW(t,z)| (3.93) 

zes.^'',\t\<T 

If T is fixed, pushing the expansion up to 2N instead of N we get easily Theorem I 1 . 2 1 using Duhamel 
formula. 

Using global estimates on aj{t,z) obtained from the transport equation (j3.76p and pushing the 
asymptotic expansion up to 2N, we get the proof of Theorem 1 1 . 41 using again Duhamel formula. 
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4 Varying phase. Proof of Theorem 11.31 

To avoid technicalities we fix the time t. It would be not difhcult to follow a time parameter t if 
necessary for application. So in this section cj) is a symplectic diffeomorphism in M.'^'^, such that cj), 
(f>~^ are Lipchitz continuous and (p G 01(20?). 

We denote z = {q,p) G K^'', (j}{z) — {Q{z), P{z)) eS.'^ xR'^ and S an action for (j), i.e a primitive 
on M.'^'^ of the closed 1-form PdQ — pdq. We consider the following phases 

4>(^'«^r) ^ + p . _ Q) _ p . (y „ 5) + 1 (e(a; _ Q) . (.^ _ g) „ r(y - g).(y - g)) (4.94) 

This class of Fourier-Integral operators with complex quadratic phase was ready analyzed in [29j . 
We want here to show how to vary the choice of the matrices 8, F for a given canonical transfor- 
mation (j> of M^'*. As in section 3, let us denote I(a, $) the operator with the Schwartz kernel 

Ka(x,y)^{2TTh)~^''/^ [ e**'*'^'"("'"'2''a(z)dz (4.95) 



where a e Co (2d), $ = ^('1'^^-^). 

Using a Fourier-Bargmann transform and the following estimate : there exist C > 0, c > such 
that for all X e R^'' we have 

\{^^,^x)\<Cexp[-^), (4.96) 

we can estimate the Fourier-Bargmann transform Ka{X,Y) of Ka and prove that X(a, $) is 
bounded in L^(K'') (see section 3, Lemma [331 and the section 5 below). 

Our goal in this section is to prove the following result which gives Theorem 11.31 as a particular 
case. 

Proposition 4.1 Let be 4 matrices in S+(d), 8,8',r,F' and a G Oo{2d). 8,8' may be z depen- 
dent such that 

3c > 0, SeWf.w > c\v\^, Vz e K^rf (4 
V7, 17I > 1,3C^, ll^jeWll < Gy, Vz e M.^'^. (4.98) 

Then there exists a semi- classical symbol a' ^ order such that we have for the Lp' 

operator norm, 

Moreover we have for the principal symbol ag the formula 

a',iz)^a,iz)^^§^ (4.100) 
° deti/2(M(0)) 

w/iere M(s) := C + DF - (^(1 - s)8 -f sO') (A + BT) 

Proof. The method is rather simple and is an extension of what we have already done for quadratic 
Hamiltonians (Proof I) except that here we have to solve transport equations in the deformation 
parameter s to get the lower order correction terms. 

Let us remark that this class of Fourier-integral operators is closed under adjointness : 

I(a, $(Q'i^))* = X(a*, $*), (4.101) 
where a*{Z) = a{(t)-^Z), Z = {Q,P), Z = (t>(z) and 

{Z- X, y) = -Sir^Z) +p-{x-q)-P-{y-Q) + 

i (r(x -q)-ix~q)- e(y - g).(y - g)) (4.102) 
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So by transitivity we can assume that F = F'. As in the quadratic Hamihonian case let us introduce 
= (l-s)e + se', ^'•"^ = $(6),,r)^ O < s < l and look for a semiclassical symbol a'-"^ = J^j ^^""^P 
such that 

d I Qi<s>'-^H^-i=^,v)ai-\z)dz = 0(ft°°), Vs e [0, 1] (4.103) 
OS 7R2d 

But we have 

(z; X, y) ~Q)-{x-Q) (4.104) 

and we have to find a C-^ family symbol a^^\ < s < 1 such that 

ll^dsa^"^ + ^ ^^^"^ - Q) • (x - Q)a('*\ $) = 0{h°°) (4.105) 
The principal term Aq = a^^^ is computed as in the quadratic case. 

Let us suppose for a moment that O, 0' are constant. Then as in the quadratic case we have 

(9, + f ap)$(") = (C" + TD^ - (A^ + TB^)es) {x - Q) (4.106) 

where A = dqQ, B — dpQ, C — dqP, D — dpP and F = is a symplectic matrix. 

We know that M{s) :— C + DT — Qs{A + BY) is invertible so we can integrate by parts as in 

section 3. and as above we can achieve the proof of Proposition 14. II 

When 8, 9' are z dependent, the integrations by part are more tricky. We have to use 

{dq + f ap)$('^) = Mt{s, z){x -Q) + N{s, z){x-Q.x^ Q) (4.107) 

where N{s^z){x, y) is a bilinear application in {x, y) E'R'^ x M.'^ into dx d matrices, with coefhcients 
in Oq in z, in s. 
Hence we have 

{x - g)e^**'''" = -^{Nrr\s, z){dq + fdp)e^'''^"' - M^)-{s, z)N{s, z){x-Q,x- Q)eT*''''"' 

(4.108) 

So we apply (|4.108p and the following lemmas to proceed like in section 3. 

Lemma 4.2 For any symbol b £ Oo{2d), for every multiindex a G N'^'' and every N > \a\/2 we 
have 

f {x-Q)"ei'^''^b{z)dz^ V ftl^l / U,0{s,z)ei'^''^d^b{z)dz + 

^<\0\<N 

h^-'^ [j^,p{s,z){x-Q)Pei-^''\p,.,dn{z)dz (4.109) 

|/3| + |7|=Ar+l,|/3|>l -'^^'^ 

where fa,p{s, z), ga.pis, z) are symbols of order 0, uniformly bounded in Oo{2d) for s G [0, 1]. 
Lemma 4.3 For every b G 00(2^) and j3 G N'' we have the crude I? estimate, uniform in s £ [0, 1], 

||I((x - Q)^b, $(^)|| = 0{h\^^^^) (4.110) 
Using these two lemmas we get the full semiclassical symbol a' a'j , where 

, , , deti/'(M(s)) ,^ 
ao(z) = ao r-r^ — ^— 4-111) 

and for j > 1, a'j is computed by induction as solution for s = 1 of the differential equation 

dsajis) = TT{M{s)M~^{s))aj{s) + bj{s), 0^(0) = a,. (4.112) 
where bj{s) depends on the afc(s), k < j — 1 □ 
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Remark 4.4 Considering the adjoint operator, it is possible to exchange the role of the matrices 
Q and T. 

If the symbol a depends smoothly on some parameter A, it is not difficult to show that a' also 
depends smoothly in X. 

Proof of Lemma 14. 2i This is done by an induction on TV such that a < N.D 
Proof of Lemma 14.31 Let us begin by giving a simple proof of (|4.96p when 8 is z dependent 
satisfying the assumptions (|4.ip . We shall prove the more general estimate, for every f3 £ N"^ there 
exist C > 0, c > such that 

< Ce-'^l'*'!', VF G M^'^ (4.113) 

Let us denote Y = {y,ri) eW^ x M.'^. By a direct estimate we get easily, 

\{x'^g'',gY)\<Ce-''^y\\ yiy,r^) eR^". (4.114) 

Using Fourier transform and Plancherel formula, we exchange y and rj and we get (|4.113p . 

Now we can follow the method of section 3 to estimate norm of operators using a Fourier- 

Bargmann transformation. 

Let be K{X, Y) the Fourier-Bargmann kernel of l{{x — QYb, <&('')). We have 

k{X, Y) = 2-^'/\^h)-'h\"\/^ f {fiZ)Anix^g''), ^Y){vx.Vz)b{z)ei'^'^^Uz (4.115) 

where Z = {Q,P) = (l>{z) and 

|(r(Z)A,(xV),^y>l - l(2:V,5-£>|. (4.116) 

So we get 

\k{X,Y)\<Ch\''\/^ f exp(-|(|y-0(z)|2 + |X-z|2)dz (4.117) 

Using that is a Lipchitz canonical transformation, we have, for Co large enough and co > small 
enough, 

\k{X,Y)\ <Coh\^\/^ e^^(-^j-{\Y ~ 4>{X)\') (4.118) 

Hence we get the proof of Lemma [4.31 using Lemma □ 
we have proved Proposition 14.11 and Theorem 11.31 

5 Semiclassical Fourier Integral Operators 

In ^23], [8] and in the recent preprint |30], the authors have considered Fourier- integral operators 
defined by the following simpler phase 

vE'(^'Q)(p;x, y) = S{y,p) + P ■ {x - Q) + ^e{x ^ Q) ■ {x ^ Q) (5.119) 

where {Q,P) = 4>{y,p), (f> is a bilipchitz canonical transformation like above, O e E_|_(c?). 
In |23| and [Sj the authors have proved semiclassical expansions for the propagator of Schrodinger 
equation for initial data with a compact support. This result is extended in [30 for the Schrodinger 
Hamiltonian — /i^A + V, to general data in with uniform norm estimates. We shall give here 
some extensions of results of [50] using the same techniques as in section 3 and 4, so we shall not 
repeat the details. 

Let us denote J^{a, ^I*"^'®) the operator whose Schwartz kernel is 

K{x,y)^{2TTh)-'^ [ e^***''"(P'^'2')a(?/,p)dp (5.120) 
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A natural question discussed in this section is to compare the Fourier-Integral operators I(a, <i)("^'®''")) 
defined with 2d "frequency variables" and J'{a,'if^'^'^^) defined with d "frequency variables". 
A Fourier integral operator in L'^(W^) is always a quantization of a canonical transformation in the 
cotangent space T*(K'^). A nice way to make clear this relationship is to use a Fourier-Bargmann 
transform (see |71[3I]). This can be easily done in the same way for Semiclassical- Fourier-Integral 
operators as we shall see now. 

Definition 5.1 A family of operators, depending on a small parameter ?i g]0, 1], [/^ : 5(R'^) ^ 
iS'(M'^) is a Semiclassical-Fourier-Integral operator of order m € R associated to the canonical 
bilipchitz transformation (j> ■ T*{R'^) T*{M.'^), if for every N' we have — U%, -f R%, where 
U%, : S{W^) -> S'{W^) and \\R%,\\ = 0{h^' ) and for every TV > there exists Cn such that 

\K''\x,Y)\ < gjvfe"^-3^/^(i+ ^^^^^^^ y^, yx,YeR^'^, he]0,l], (5.121) 

where K'^{X,Y) is the Schwartz kernel of J-bU^iJ-^. 

Remark 5.2 1. In this definition, which coincides with a definition given in J31f for h= 1, a 
semiclassical-Fourier-Integral operator has, up to a negligible operator in h, a kernel living in 
a neighborhood of the graph of a canonical transformation cj). But this definition says nothing 
concerning asymptotic expansion ofK^{X,Y) in a neighborhood of the graph of (j) when H is 
small. So this definition is certainly too permissive. But for h fixed it is suitable as proven 
in fSIF- 

2. Using Carleman-Schur estimate, a semiclassical F.I.O of order is uniformly bounded in 
L^(]R''). This is a straightforward consequence of the definition. This class of semiclassical 
F.I.O of order is clearly closed by composition. 

/p\ 

3. In Definition 15.-?! it is equivalent to use any Fourier-Bargmann transformation J-j^ , F g 
S+(d). 

4- There are other definitions of semiclassical F.I.O using Lagrangian analysis and real phase 
functions. For this point of view see for example fl^. 

5. Fourier- Integral operators with complex phase were used to study propagation of singularities 
of P.D.E. Many papers and books have been published on this subject, among them let us 
point out \26, 32]. 

Now we shall see that the operators already considered in this paper are semiclassical-Fourier- 
Integral operators. 

Proposition 5.3 Let be amplitudes a = a{x,z), a £ Oo(3d) and u — u{x,y,p), u G 00(3^) 
and 0,F G I]-|_(d), may depend in z or {y,p), such that il.l6\) . Ijl.lT^ are satisfied. Then 
X(a, <I'^'^'®''"-') and J^iu,'^'^'^) are semiclassical-Fourier-Integral operators of order 0. 

Proof. Concerning X(a, $*^'^'®''"^), we get the result following subsection 3.2, estimate (|3.92p . 
The proof for J^{u, ^I*"^'®) is almost the same. For simplicity we assume 6 constant. For 8 
depending in (jj,p) we could proceed as in section 4. 
Let us denote X = {x,^), Y — {y, f}). We want to estimate 

k{X,Y) [2TTh)-'^ 1 1 1 ei^u{x,y,p)dpdxdy (5.122) 

where 

$ = Siy,p) + P-ix^Q) + ^{x ~ Q) ■ {x -Q) + 

-{i~y)-{S:-y)+^-{S:-y) + -{y-x)-{y-x) + fi-{y-x) (5.123) 
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Let us remark that we have : F'^^ = ( i^'^ ^ ~ \c Di 

symplectic, we know that D'^ — B'^Q is invcrtible. Hence we have 

dy^ = {C^ - A''&){x-Q) + {i-p) + t{S:-y) (5.124) 

= {D'' - B''Q){x~Q) (5.125) 

= e(a;-Q) + (P-r7)+i(y-a;). (5.126) 

So we get the necessary estimates on K by integrations by parts using 

dy^-{-A^e + c^)iD^ -B^e)-^dp-9 ^{i^p) + i{i-y) (5.127) 

da:^ -0(0'^ - B'^ey^dp'i' ={P -f]) + i{y-x) (5.128) 

□ 

The foUowing resuh is a shght generahzation of pSl [HI [30] ■ 
Theorem 5.4 Under the assumptions of Theorem \1.2\ and U.16]) . |J. j7p , we /laiJe 

K{t;x,y)^{2TTh)~'^ [ e*'''^''' ^^'^•P''''^u{h;t,y,p)dp (5.129) 



where u{h;t,y,p) = Uj{t;y,p)h-' has the same meaning as in Theoren il.Si 

0<j< + oo 

/n particular 

uoit,y,p) ^ det^/^(D -OB)) (5.130) 



Sketch of proof. These resuh can be proved foUowing the same strategy as for proving Theorem 

We first prove the Theorem for some 9 (0 = if), following the proof of Theorem 1 1.21 then we can 
get the Theorem for any 8 by the variation argument as in the proof of Theorem 11.31 estimate 
for operator norm of Fourier- Integral operators is used to control the remainder terms. □ 

Remark 5.5 It is not difficult to adapt the proof of Theorem \1.4\ concerning an Ehrenfest time 
estimate to the setting of Theorem \5.4\ 
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